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The creation of ultracold molecules is currently limited to diatomic species. In this letter we 
present a theoretical description of the photoassociation of ultracold atoms and molecules to create 
ultracold excited triatomic molecules, thus being a novel example of light-assisted ultracold chem¬ 
ical reaction. The calculation of the photoassociation rate of ultracold Cs atoms with ultracold 
Cs 2 molecules in their rovibrational ground state is reported, based on the solution of the quantum 
dynamics involving the atom-molecule long-range interactions, and assuming a model potential for 
the short-range physics. The rate for the formation of excited Css molecules is predicted to be com¬ 
parable with currently observed atom-atom photoassociation rates. We formulate an experimental 
proposal to observe this process relying on the available techniques of optical lattices and standard 
photoassociation spectroscopy. 


Ultracold dilute gases at temperatures much lower 
than 1 mK offer new opportunities for the study of ele¬ 
mentary reactive processes between atoms and molecules. 
As emphasized in seminal review articles [Ei this so- 
called ultracold chemistry is freed from averaging on large 
velocity distributions. Thus in this regime the quantum 
nature of the processes can be accessed, like the presence 
of resonances, or the sensitivity to long-range interac¬ 
tions which can induce anisotropic arrangements prior to 
the reaction. One of the most advanced experiments in 
this direction is undoubtedly under progress in JILA at 
Boulder, where quantum threshold collisions induced by 
tunneling through a centrifugal barrier between fermionic 
ultracold KRb molecules have been detected ii- Many 
experiments demonstrated that conditions suitable for in¬ 
ducing chemical reactions between ultracold ground state 
alkali-metal atoms and molecules tightly bound in their 
electronic ground state are reachable. Ultracold collisions 
between Cs atoms and Cs 2 [11,0 or LiCs 7| molecules, of 
Rb and Cs atoms with RbCs molecules [8| , or of K and 
Rb atoms with KRb molecules 0 have been detected 
through atom trap losses. Three-body recombination in 
a Rb quantum gas has been characterized 0 , and fea¬ 
tures associated to universal A^-body (up to = 5) reso¬ 
nances induced by collisions between atoms and weakly- 
bound molecules have been observed in quantum degen¬ 
erate gases 

However the deep understanding of the collisional dy¬ 
namics is still to come as none of these achievements were 
able to characterize the nature and the state distribution 
of the final products EE!. The heavy mass of alkali- 
metal species leads to a huge density of resonant states 
related to the numerous rovibrational levels accessible 
during an atom/molecule or a molecule/molecule colli¬ 
sion, preventing them to be individually characterized. 
Models connecting the standard treatment of long-range 
interactions between particles to a statistical description 
of the resonances mostly resulting from short-range cou¬ 
plings have been developed |18l42l|. The main quantita¬ 


tive uncertainty of such models arises from the estimation 
of the actual amount of resonant states which are indeed 
active during the collision. 

In this letter, we propose to use the well-known pho¬ 
toassociation (PA) process to get more insight into such 
collisions between ultracold atoms and molecules. PA of 
an ultracold atom pair is a laser-assisted collision creat¬ 
ing an electronically-excited diatomic molecule often in 
a weakly-bound energy level 22, [2^. Under favorable 
circumstances such a short-lived molecule can decay by 
spontaneous emission (SE) into a rovibrational level of a 
stable electronic state. The PA-j-SE process played a cru¬ 
cial role in creating samples of stable ultracold molecules 
23-22 1- This is the simplest example of the formation 


of a chemical bond under ultracold conditions. The PA 
step is mostly controlled by long-range atom-atom inter¬ 
actions |28I43QL while the SE step involves short-range 
interaction |24 EE. The knowledge of atom-atom in¬ 
teractions is thus required to fully describe the PA-j-SE 
process, which is possible in most cases. 


Here we discuss the PA of ultracold alkali-metal atoms 
and diatomic molecules to create ultracold trimers. We 
evaluate the rate for PA of ultracold ground state Cs 
atoms and ultracold ground-state Cs 2 molecules in their 
lowest rovibrational level, revealing a magnitude similar 
to the rates observed for atom-atom PA. Our model only 
involves the long-range interactions between the parti¬ 
cles, thus neglecting the influence of the short-range in¬ 
teractions responsible for a possible high density of reso¬ 
nant states. Thus the problem is formally similar to the 
PA of a pair of atoms of different species. We propose 
an experimental way to observe the excited Csg molecules 
using an nearly-degenerate quantum gas of Cs atoms and 
molecules trapped in an optical lattice [HI • By compari¬ 
son to the present simplified model, such an observation 
would shed light on the link between short-range and 
long-range physics in ultracold processes. 

We consider the interaction between two charge dis¬ 
tributions A and B with total angular momentum Ja 









and Jb and separated by a distance R along an axis 
Z joining their center-of-mass and oriented from the 
molecule toward the atom. It can be split according to 
two distance ranges: (i) the long-range domain where 
interactions are determined by the individual properties 
of A and B and can be accurately calculated, and (ii) 
the short-range domain where complex chemical inter¬ 
actions take place. Both domains connect around the 
LeRoy radius Rlr accounting for the size of each dis¬ 
tribution (ssj. For the rest of the paper we consider a 
Cs( 6 ^Si/ 2 ) atom and a Cs 2 molecule in the lowest vi¬ 
brational {v = 0 ) and rotational level A" = 0 of its 
electronic ground state as the initial state \i) of 

PA. We have R^^^ = 42ao (ao = 0.052917721092 nm). 
The laser frequency is assumed to be slightly detuned 
by 6pA from the D 1 or D 2 Cs line, so that PA popu¬ 
lates energy levels of the CS 3 complex located below the 
Cs( 6 ^P 3 / 2 )+Cs 2 (X^I;+,i; = 0 ,A = 0 ) dissociation limit. 
For these final states we have Rpp = 46ao. We ignore 
the Cs( 6 ^Pj) hyperfine structure for simplicity, and we 
will consider experimental conditions where it can be a 
secure approximation. 


Following our previous work, the long-range excited 
atom-ground-state molecule interaction is treated within 
the second-order degenerate perturbation theory 3^ 35[ . 
In brief, the Hamiltonian of the system is written as 
H = Ao+W, where Hq refers to the energy of the individ¬ 
ual particles at infinity. The atomic states are labeled as 
|/ = 1, s = 1/2, j, cj), where j = 1/2, 3/2 is the total angu¬ 
lar momentum quantum number and u its projection on 
the Z-axis. The molecular states |X^S+,'i; = 0, A, ruAr) 
involves the projection m at of A onto Z. The mutual ro¬ 
tation of the atom and the molecule is not introduced yet, 
so that the total projection quantum number Q = mAr+cj 
characterizes the trimer states which will be referred to 
as |jf, cj. A, mAr; f2;p), where p = (—1)^ is the parity. 
The W(R) operator includes the first-order qnadrupole- 
quadrupole term (Vqq{R) oc 1/R^, if A 7 ^ 0) [s^, and the 
second order dipole-dipole interaction {vj^\R) oc 1/R^) 


The Hamiltonian written in the |j, cj. A, ruAr; ba¬ 
sis is diagonalized, yielding the adiabatic potential en¬ 
ergy curves (PECs) for large atom-molecule distances 
displayed in Fig.[Tl for j = 1 / 2 , 3/2, = 1 / 2 , 3/2, and in¬ 

cluding rotational levels up to A = 5. The complex struc¬ 
ture of the PECs is due to the coupling between the pure 
electrostatic interaction and the rotation of the molecule. 
In particular, these PECs cannot be expressed anymore 
as a pure R~^ expansion. Eor instance the lowest curve 
labeled with A = 0 at large distance for j = 3/2, (pan¬ 
els c and e of Eig. [T]) is more attractive - and thus more 
favorable for long-range PA- due to the coupling with the 
A = 2 state (see Eig. 2 in Ref. Hi). We also computed 
the PECs including only A = 0 , 2 levels in the Hamilto¬ 
nian, showing that upper rotational levels only slightly 



FIG. 1. (Color online) Long-range potential energy curves 
between a ground-state Cs 2 (X^S^,i’ = 0,A) molecule and an 
excited atom Cs(6^Pj) as a function of the scattering coordi¬ 
nate R. (a),(b): j = 1/2, Q = 1/2; (c), (d): j = 3/2, Q = 1/2; 
(e), (f): j = 3/2, Q = 3/2. The origin of energies is taken 
at the Cs 2 (X,z/cZ = 0,A=0)-hCs(6^Pj=i/2) limit. Closed cir¬ 
cles in (a) and (c) correspond to PECs computed with only 
A = 0, 2 levels in the Hamiltonian. 


affect the lowest PECs for A = 0 . At R ^ Rlr the 
PECs are attractive by at most 1 cm“^, so that we limit 
our study to PA detunings smaller than 1 cm“^ below 
the Cs D 2 line. Note that the PECs correlated to the 
Cs( 6 ^Pi/ 2 )+Cs 2 (X^E+,i; = 0,A) are obviously simpler 
than in the j = 3/2 case, which could yield preferable 
conditions for a first experimental implementation. 

Eor R < Rlr, each diagonal term of the above Hamil¬ 
tonian is smoothly matched to a Lennard-Jones poten¬ 
tial Vlj = Dlj{Clj/R^){{Clj/R^) - 1). The Clj 
factor is chosen equal to the computed Cq coefficient 
[ 3 ^ . while the depth Dlj is taken large enough to 
ensure that the short-range part of the wavefunctions 
have a negligible contribution to the ultracold dynam¬ 
ics. Eor the initial PA state we have Cq = 12101 a.u. 
and Dl = 10867.2 cm“^. Eor the final PA states of 
interest here, namely with j = 3/2, A = 0,2 and 
\Q\ = 1/2, the parameters are Cq{N = 0) = 60311 a.u., 
DUN = 0) = 1500 cm-\ and C/(A = 2) = 25311 a.u., 
Df{N = 2) = 500 cm-h 



















































The radial wavefunctions f{Cl,v')) associated to the 
eigenvalue Ef{v') of the coupled states described by the 
above PECs are calculated with the mapped Fourier grid 
Hamiltonian (MFGH) method [38| suitable for calcula¬ 
tions at large distances as the grid step is mapped onto 
the local kinetic energy of the channels. We used a grid 
of Ng = 4145 points between 4ao and TOOOao with the 
MFGH scaling parameter of (3 = 0.11 [HI to ensure a 
proper convergence of the Ef{v') values. 

We assume that PA takes place during a s-wave atom- 
molecule collision at an energy Ei = ksT^ where ks 
is the Boltzman constant and T the related tempera¬ 
ture. The associated continuum radial wavefunction 
is computed by solving the Schrodinger equation for the 
entrance channel with the Numerov method. In the per¬ 
turbative regime, if the PA laser with intensity /pa, de¬ 
tuning ^PA and polarization epA, hits a level v' of the Csg 
complex, the PA rate is expressed as [HI 
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( 1 ) 

where Ap = ^ fi is the atom-molecule 

reduced mass, rimoi stands for the Cs 2 density, and 
^ifi^PA) is an angular factor depending on ^pa- Its value 
strongly depends on the initial state of the particles, and 
we took Aif = 1/2 as a typical example. As the detuning 
is assumed to be small, the \i) \f) transition dipole 

moment is taken as the one for the 6^5'i/2 ^ 
transition, leading to the corresponding atomic Rabi fre- 
quency flif and saturation intensity Iq = 1.1 mW/cm^ 
[28[. For the purpose of comparison among various ex¬ 
perimental implementations, it is useful to define the nor¬ 
malized PA rate Kpa{T) = RpA{T)/nmoi/4^PA (conve¬ 
niently expressed in cm^) where 0 pa = IpA^PA/hc is the 
PA laser photon flux at the wavelength Apa [13,01 • The 
results for Kpa{T) are displayed in Fig. [2^ for two differ¬ 
ent temperatures: at T = 20 fiK typical for a magneto¬ 
optical trap (MOT), and at T = 500 nK representative 
of a gas close to the quantum degeneracy (QD) regime. 
Like in the case of heteronuclear diatomic molecules 0, 
the oscillating pattern reflects the strong variations of the 
spatial overlap between the initial and final radial wave- 
functions with the detuning. As expected from Eq. o, 
PA close to the QD regime and for small detunings ex¬ 
hibits a large r ate, as already demonstrated in the exper¬ 
iments on Li 2 411 and Na 2 0. We mention that the 
coupling of the long-range multipolar interactions with 
the rotational energy quoted above also contributes to 
enhance the PA rate as the entrance channel A" = 0 is 
more attractive. This first principal result is further il¬ 
lustrated in Fig. [213 showing that the predicted PA rate 
in a MOT environment is comparable to the rates ob¬ 
served with various heteronuclear diatomic species. We 
recall that the PA rate for Cs 2 at 140 /iK was larger 
than most of the PA rates observed under MOT condi- 


FIG. 2 . (Color online) (a) Normalized PA rate ATpa as a 
function of the detuning Spa of the PA laser with respect 
to Cs 2 (X^S^, 7 ; = 0,A = 0) -h Cs( 6 ^P 3 / 2 ). The points 
correspond to vibrational levels v' of Css for |Q| = 1 / 2 . 
The rates are calculated at T = 20 /jK (open red circles) 
and T = 500 nK (closed blue circles). (b) Comparison 
of computed Kpa values at 140 /xK and 500 nK for two 
CsWibrational levels v' at 0.770 and 1.008 cm“^ with the 
theoretical value in Cs 2 [ 1 ^, and several experimental rates 
for LiK in, LiRb 0, RbCs 0, Naa 0, and Lb 0. 

Although the results look very promising, we must 
discuss whether the conditions of density and tempera¬ 
ture for atoms and molecules needed for the formation 
of Csg can be satisfied with the available experimen¬ 
tal techniques. A density rimoi of Cs2(X^F+,z/=0,A=0) 
molecules of about ~ 10^^ to 10^^ cm“^ and temperature 
T ^ 500 nK, can be obtained departing from an ultra¬ 
cold sample of ground-state Cs atoms [32| with similar 
density. These are favorable conditions to initiate a PA 
experiment where products would be indirectly detected 
through the loss of molecules. 

As discussed above, the direct detection of products 
is also a central objective to better understand the full 
dynamics of the process, which could benefit from the 
presence of an optical lattice (OL). The scheme of our ex¬ 
perimental proposal is presented in Fig. [3l based on the 
conditions of Ref. [HI • The motional degrees of freedom 
of the atoms are controlled by an OL with wavelength 
Aol=1064.5 nm. The OL intensity is tuned in order 
to induce a Mott-insulator (MI) state with preferentially 
two atoms per lattice site (see e.g. Ref. [i^). The magne¬ 
toassociation technique is then applied, creating ground- 
state Cs 2 molecules in a weakly-bound vibrational level. 
This population is transferred to the Cs 2 rovibrational 
ground state by means of stimulated Raman adiabatic 
passage (STIRAP) technique which yields a molecular 
MI state to a good approximation [32| (Fig. [3^). 
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During this sequence a fraction of ultracold atoms re¬ 
mains unpaired. We propose to tune the OL intensity to 
drive the transition from a MI phase to a superfluid (SF) 
phase for atoms, keeping the molecules in the MI state 
(Fig. [Sb). Assuming a deep optical lattice with spacing 
between two successive sites d = XolI‘^^ the MI-to-SF 
phase transition for the species /3 (/3 = Cs or CS 2 ) oc¬ 
curs when the lattice depth Vq oc cx^Iql (where /ql is 
the corresponding intensity and the dynamic dipole 
polarizability of the species at Xql) satisfies the critical 
value given by 



{Vo% \ 
E^r ) 




( 2 ) 


where is the recoil energy of the considered species 
and is the corresponding scattering length. The crit¬ 
ical value /J^)c of the ratio between the on-site in¬ 
teraction and the tunneling energy depends on 

the number of (3 particles per OL sites n. We will as¬ 
sume n=l hereafter, leading to (JJj = 29.36 [3^ in 
the case of a simple cubic lattice. At Aql = 1064.5 nm, 
the recoil energies are such that jkB = 64 nK and 
E^^^lks = 32 nK. Equation (|2]) shows that the criti¬ 
cal value depends on the scattering length of the two 
species (taken as positive quantities as a MI phase can 
be achieved) but only through a logarithmic way. While 
the variations of are well established [T3, l48| , those 
of totally unknown. A reasonable assumption is 

to consider that the experiment is performed under con¬ 
ditions where and a^®^ are not taking vanishing or 
infinite values, such that d/a^ ^ 1. This leads to the 
estimation of the ratio of the critical OL intensities for 
each species 

{iSDc^ 

which must be smaller than unity to ensure that the MI- 
SF transition is carried out on the atoms without re¬ 
leasing the MI phase of the molecules. Assuming that 
^ 2(a^® [ 43 , the ratio indeed amounts to about 
1/4. This estimate is further confirmed in Fig.jSJi, show¬ 
ing that (^OL )c/(^ol)c is indeed smaller than unity for 
a broad range of scattering lengths (Fig. [3t) for both 
species. 

The PA laser is not expected to significantly disturb 
such an hybrid system. Indeed, for Apa ~ 852 nm (or 
^ 11737 cm“^, close to the D2 Cs line) the Cs 2 ground 
state molecule will only weakly absorb such a light due to 
unfavorable Franck-Condon factors with excited molecu¬ 
lar states, hence the related scattering force will be neg¬ 
ligible (see Fig. 3 of Ref. 0). 

Finally, assuming that the photoassociated Csg com¬ 
plex has a radiative lifetime similar to the one of 
electronically-excited Cs or Cs 2 species, the products 


FIG. 3. (Color online) Scheme of the experimental proposal. 

(a) Preparation of a MI state for both atoms and molecules. 

(b) The PA laser is applied to a mixed molecular MI - atomic 
SF state, (c) Excited CS 3 molecules are created, and trapped 
after radiative decay into stable Cs 3 trimers. (d) Contour 
plot for the ratio of critical OL intensities for Cs and Cs 2 
(Eq. dS])) as a function of their respective scattering lengths. 
The dashed line indicates where this ratio equals unity. 


of the decay could be discriminated in the experiment 
(Eig.[3t). On one hand, the decay of CS 3 back into three 
Cs atoms or a Cs/Cs 2 pair with a kinetic energy release 
will result in a loss of the fragments from the OL. On 
the other hand, if stable Css molecules are created, they 
will most probably be trapped in the OL, as their recoil 
energy is even smaller than the Cs 2 one as their larger 
dipole polarizability induces a deeper trap. A branching 
ratio between the two decay channels could thus be in¬ 
ferred from the benefit of the ” ‘half-collision” ’ character 
of the process leading to bound products. 

The observation of ultracold trimers in a PA experi¬ 
ment will shed light on the role of short-range interac¬ 
tions in ultracold collisions in between two extreme sit¬ 
uations: either PA results into well identified resonant 
features assigned to Csg levels determined by pure long- 
range interactions (as in the present model), or PA will 
reveal a quasicontinuum of levels induced by the strong 
coupling of such long-range states with numerous reso¬ 
nances of the trimer (as in the approach of Refs [30. 
The presented approach can be generalized to stable ul¬ 
tracold Cs 2 molecules in their lowest state, but this 
case could be quite challenging as the PA wavelength will 
reach the region of numerous Cs 2 excited levels (see Eig. 6 
in Ref. 0). Other species can also be investigated, as 
for instance the association of Li atoms and K 2 ground 
state molecules where the K 2 rotational energy is similar 
to the Li* spin-orbit splitting, inducing complex patterns 
in the long-range PECs (ssj. 
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